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Abstract 

Within a field theoretical formalism suited to treat inhomogeneous hot 
quantum systems, we derive the two-particle correlation function for particles 
having a spectral width 7 in the region of their emission. We find, that this 
correlation function measures the radius Rq of the thermal source only in case 
7^0 > 1. 
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In the field of relativistic heavy-ion collisions the analysis of Bose-Einstein cor- 
relations has attracted much attention [I], 0, ||J]. The general hope is to extract 
information about the size of a source radiating pions by studying their two-particle 
correlations. These correlations are typical quantum effects, hence quantum field 
theory is the proper framework to address the problem theoretically. Two types 
of sources are physically interesting: Classical currents and thermal sources - the 
physical reality being a mixture of these cases. 

The problem of a free quantum field radiating from a classical current is exactly 
solvable pp.438], also at finite temperature |J. Conversely, it is quite difficult 
to find a consistent theoretical description for bosons radiating thermally from a 
local "hot spot". One reason is, that such a physical situation corresponds to a 
non-equilibrium state. Another reason for this difficulty is, that in thermal states a 
breakdown of perturbation theory may occur when it is expressed in terms of stable 
particles having zero width @. 

With the present paper we attempt to study these two aspects, i.e., we address 
the questions: 1. How are Bose-Einstein correlations affected by gradients in the 
temperature distribution, and 2. How are they affected by a finite lifetime of the 
bosons. 

In our model, we assume an entirely thermal radiation of pions from a "hot 
spot", which acquire a nonzero spectral width in the region of their generation, e.g. 
by coupling to A33-resonances. All thermally generated pions are propagating out 
of the hot spot. 

This picture is grossly simplified as compared to the physical reality: At least 
a substantial fraction of pions generated in relativistic heavy-ion collisions is not 
produced thermally, but arises from the decay of the A 33 -resonance outside the 
"hot spot". However, on one hand it is legitimate to study only a partial aspect 
of the full problem. On the other hand, the formalism introduced here is general 
enough to be extended to the resonances as well. 

Quantum field theory for non-equilibrium states comes in two flavors: The 
Schwinger-Keldysh method || and thermo field dynamics (TFD) ||. For the pur- 
pose of the present paper, we prefer the latter method: Apart from its technical el- 
egance, the problem of an inhomogeneous temperature distribution has been solved 
explicitly in TFD, up to first order in the temperature gradients |TI], FT]. This 
solution includes a nontrivial spectral function of the quantum field under consid- 
eration: It employs a perturbative expansion in terms of generalized free fields with 
continuous mass spectrum M. 



described by a statistical operator (or density matrix) W, and the measurement of 
an observable £ will yield an average that is calculated as the trace of £W over the 
Hilbert space of the system. In thermo field dynamics (TFD), this is simplified to 
the calculation of a matrix element 

l ' V ((H|w--)) w 

and the Hilbert space is doubled ||. The thermal pion field is described by two 
scalar field operators 4> x , <fi x and their adjoints <f>*, (f>*, with canonical commutation 
relations. The field <p x is expanded into momentum eigenmodes: a^(t) creates a 
pion with momentum k and charge I = ±1, a second set of operators (commuting 
with a) \ a) exists for the tildean field 6„ H, ITT . 



These operators do not excite stable on-shell pions. Rather, they are obtained 
as an integral over more general operators £, £ with a continuous energy parameter 
E: 



"®) = JdE fa* q A- /e (£,(q + k)/e) (B-(£, q ,k) 
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The principles of this expansion have been derived in ref. 0, A(£, k) is a positive 
weight function. For equilibrium states, this function is the spectral function of the 
field <p x . For non-equilibrium systems, the existence of a spectral decomposition 
cannot be guaranteed 0. We may expect however, that close to equilibrium the 
field properties do not change very much. Thus, with this formalism we study a 
quantum system under the influence of small gradients in the temperature, with 
local spectral function A(£, k). Corrections to such a picture only occur in second 
order of temperature gradients [|TI], |TT|j . 

A thorough discussion of the 2x2 Bogoliubov matrices was carried out in ref. 



12]. For the purpose of the present paper, we simply state their explicit form as 



B,(E, q ,k) = \ «"(«-*) + ".<*,,,*)) -*<*,«,*), (3) 

1 -6 3 (q-k) 6 3 (q-k) 1 

where N(E, q, k) is the Fourier transform of a space-local Bose-Einstein distribution 
function 

N t (E,q,k) = fd^ze-^- h >n l {E,{q + k)/2,z) 



The ^-operators have commutation relations 

~Zeu, d'm] = Sw 5{E - E') 6\k - k') . (5) 

Similar relations hold for the £ operators, all other commutators vanish, see 0. 
They act on the "left" and "right" statistical state according to 

Z Ek i\\W)) = 0, iBulW} = 0, ((1||£L = 0, ([1||3h = VE,k,l = ±l. (6) 

With these rules, all bilinear expectation values can be calculated exactly. Higher 
correlation functions have a perturbative expansion in the spectral function. 

Of these, we are interested in the two-particle correlation function, which is the 
probability to find in the system a pair of pions with momenta p and q: 



Cii'{p, q) = -j— ^—j— r = 1 + Oil' , . j-t , (7) 



For simplicity, we abbreviate the mean momentum of this pair by Q = (q + p)/2. 
Using the above rules of thermo field dynamics, the functions and JF g are cal- 
culated as 

oo 

^oo(p) = JdEjd 3 zAi(£,p)\ t (S,p,z) 



oo 

F E (p,q) = JdE J d 3 z (A t (£,p)A t (£,Q))^ e i{p ~ q)z ni(E , Q , z) . (8) 
o 

Before we use the above expression to obtain numerical results, we perform an expan- 
sion of jF e around the mean momentum Q. This is consistent with the restriction, 
that the spectral function acquires corrections in second order of the gradients: 



Fe{p,q) = F' e (p,q) 

oo 

+ JdE Jd 3 ze^ z (^V Q A x (£,Q)V z \ l (£,Q,z) S j +C(Vf\) 
o 

oo 

r e (p,q) = JdE Jd*ze^-v>A x {£iQ)\ x {£,Q,z). (9) 
o 

The first part of this expansion is, apart from the folding with the spectral function, 
also obtained in other calculations of the correlator [0, |3|, || . This standard expression 
for the correlation function therefore is 



and our result for cu> (p, q) differs from cw (p, q) in first order of gradients in the 
distribution function. For the purpose of generalizing our result it is worthwhile to 
note that the gradient term in (^|) is just one half of the Poisson bracket of A and 
n [0. 

We have calculated these correlation functions with a simple parameterization 
of the pion spectral function, 

= ^r w^h^ (») 

v 

where D, p = a/ m 2 + p 2 + 7 2 and m w = 140 MeV. This parameterization has been 



motivated and related to a more serious field theoretical approach in ref. [13|]. To 
gain information about the maximal influence exerted by the occurrence of a nonzero 
spectral width, we studied only the case of an energy and momentum independent 
7 equal for both charges. The temperature distribution was taken as radially sym- 
metric gaussian 

T(z)=T(r) = T e*p(-^ , (12) 

with chemical potential \i — and Rq = 5 fm. The local equilibrium pion distribution 
for a given momentum k is obtained by folding n with the spectral function. Hence, 
the mean radius of this particle distribution function acquires a 7-dependence. We 
define the rms radius orthogonal to the direction of Q as 

00 00 

R n = Ij = Jdrr 1 JdEA(£, k) (e £/r ^ - oof °° . (13) 



Note, that R n is not the 3-dimensional rms radius of the distribution function (which 
would be l^j I2). Rather, R n is half the product of angular diameter and distance 
between detector and source. A constant temperature over a sphere of radius Ro 
would yield an R n = R /\/3, while its 3-D rms radius is -R a/3/5. 

The correlation functions, c«/ (p, q) calculated according to (0) and cw (p, q) as 
defined in eqn. (|T0"D then may be fitted by a gaussian form, i.e., 

c w (p, q) « 1 + exp {-R\p - qf) (14) 

and similarly for cw (p, q) with parameter R. We assume for our conclusion, that 
cw(p,q) is the correlation function measured experimentally. 

In the figure, we have plotted the two fit parameters R, R and R n as function 
of 7. The principal result of the calculation is, that for small values of 7 the 




Figure 1: Correlation radius of a pion "hot spot" with temperature T =100 MeV 
Long dashed line: Gaussian fit to Q//(p, q) according to (0) 
Short dashed line: Gaussian fit to cu>(p,q) according to (p~0]) . 
Continuous line: mean radius of the thermal distribution function, eqn. (|13|) . 



Consequently, the measured correlation radius R is always larger than R as expected 
from the function cu/(p,q). The deviation is such that for small enough 7, R « 
R+l/j. For larger 7, the small differences between R, R and R n may be attributed 
to our use of a gaussian temperature distribution: n(T(r)) is not stricly gaussian, 
only in the (unphysical) limit 7 — > 00 one reaches R = R = R n = R0/V2. 

Before we interpret this result, we have to admit that our calculation is very 
crude: Neglecting energy and momentum dependence of the 7 in the spectral func- 
tion can be a first step only Also, in a more realistic calculation the partially coher- 
ent production of pions would have to be taken into account, forcing cw{p,p) < 2. 

Within this limitation however we feel safe when stating the following answer to 
the questions raised in the introduction: A finite lifetime or nonzero spectral width 
7 of the bosons is essential, if one wants to infer the thermal source radius Rq from 
correlation measurements. To be more precise, only for 7.R0 > 1 the correlation 
function measures the mean radius of the particle distribution function. 

the 



equilibration rate of a distribution function is, to lowest order, given by the spectral 
width of the particle |T2|]. Consequently, a very small 7 corresponds to a system 
that does not equilibrate - hence the correlation function approaches the "quantum 
limit" Cu>(p, q) — > 1 + Sai5 pg , and the correlation radius obtained by a gaussian fit 
becomes infinite. 

Also, for a given energy, I/7 is a measure for the spatial size of the pion "wave 
packet", which must be smaller than the object to be resolved. In stars emitting 
photons having only their thermal width 7 « aT « 0.5/137 eV, the condition 
7-Ro 3> 1 is always satisfied. For the typical thermal pion sources occurring in 
relativistic heavy-ion collisions however, this criterion can be violated: Pions of 100 
MeV momentum in nuclear matter have an effective 7 of only a few MeV ]13| . 

We conclude, that in the Hanbury Brown-Twiss analysis of relativistic heavy-ion 
collisions one may measure correlation radii, i.e., by gaussian fits to the correla- 
tion function, which are much larger than the actual radius of the thermal source. 
Most certainly it is not possible to infer the size of a thermal source by correlation 
measurements, when the spectral function of the measured bosons in the region of 
their generation is unknown. As a rule of thumb we suggest, that pion correlation 
functions should be measured for momenta of \Q\ ~ 350 MeV, since then the mixing 
of pions with A 33 -resonance/nucleon-hole excitations is largest. 

As a final note we emphasize again, that in common calculations of the correla- 
tion function one has to introduce ad-hoc random phases between several classical 
sources and then obtains the "standard" correlator cw(p,q). Instead, we relied on 
a proper field theoretical treatment correct up to first order in gradients of the tem- 
perature. We found, that the non-equilibrium character of the system must be taken 
serious when calculating the correlation function. 
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